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We study the anisotropies of the Galactic confusion noise background and its effects on LISA data 
analysis. LISA has two data streams of the gravitational waves signals relevant for low frequency 
regime. Due to the anisotropies of the background, the matrix for their confusion noises has off- 
diagonal components and depends strongly on the orientation of the detector plane. We find that 
the sky-averaged confusion noise level S{f) could change by a factor of 2 in three months, and 
would be minimum when the orbital position of LISA is either around the spring or autumn equinox. 
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^ ; I. INTRODUCTION 

o : 

fSJ , The Laser Interferometer Space Antenna (LISA) is planned to be launched around 2011 and is expected to establish 

5^ a new window in the low frequency gravitational wave astronomy from O.lmHz to lOOmHz Its main astrophysical 
targets are Galactic binaries and cosmological massive black holes (MBHs) that merge with other MBHs or capture 

^> ■ compact objects. As for the Galactic binaries, there would be a lot of sources in the LISA band. For example we 

\ will be able to resolve several thousand close white dwarf binaries at / > 3mHz 0, ■ At lower frequency regime 

. / 3mIIz they are highly overlapped in the frequency bins, and it would be difficult to resolve them individually. 

^ As a result, they form a confusion noise background whose magnitude could be larger than the detector noise at 

►-j^ O.lmHz < / < 3mHz The coalescence frequency /c of a MBH system is given by its redshifted total mass as 

Tjj- ' /c ~ 2(M^/2 X 10'^MQ)-imHz. The mass function of MBHs is highly uncertain at the lower end ~ IO^Mq O, and 

I we might have to search treasurable signals from cosmological MBHs in the Galactic confusion noise background. 

' The spatial distribution of the Galactic binaries would trace the Galactic structure well, and the confusion noise 

\ background would be strongly anisotropic. The background itself can be regarded as a signal that would provide 

^-r > some information on our Galaxy |al2lla- But we should notice that there exists a more straightforward approach to 
probe the Galactic structure using thousands of resolved binaries whose three dimensional positions are estimated in 
some error boxes In this paper we will focus on the role of the background as a noise, and study the effects of its 

^ • anisotropies on the signal analysis of LISA. We do not pay attention to the normalization of the background that has 

1^' been studied by other papers using an isotropic approximation of the source distribution f3|. 

^ • This paper is organized as follows. In Section 2 we study the response of interferometers to the anisotropic 

rr| gravitational wave background with using the long wave approximation, and formulate the noise matrix for the 

^ . Galactic confusion background. In Section 3 we numerically evaluate the bakground noise as a function of the 

' orientation of the detector, and estimate the annual modulation of the noise level for LISA. In Section 4 we discuss 

■ the anisotropies of the signal to noise ratio (SNR) of sources at a fixed distance. Section 5 is devoted to a brief 

■ summary of this paper. In the Appendix the confusion noise matrix is analyzed in a different manner from Section 2. 



II. FORMULATION 



From data streams of LISA, we can make three meaningful modes (^, E, T) such that the laser frequency noise can 
be reduced enough and their detector noises do not correlate (see |9| and references therein, also [13 )• At the low 
frequency region (/ ^ lOmHz: determined by the arm-length of LISA, 5 x lO^km) relevant for the present analysis of 
the Galactic binary confusion noise, the T-mode has much worse sensitivity than the rest. Therefore we only discuss 
A and E modes here. These two data streams di {I = A, E) are made by the detector noise nf and the response to 
gravitational waves hj 

di = hi + nf. (2.1) 
We define the detector noise spectrum Sfj as follows 0, 

{(nfifyrnfi/)) = hif - f')Srj{f). (2.2) 
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The frequency dependence is irrelevant for our analysis, and, hereafter, we omit its explicit dependence for notational 
simplicity. We have the following relation for the detector noise spectrum matrix 

^AA = ^EE^ ^AE ~ 0- (2-3) 

The three spacecrafts of LISA form a nearly equilateral triangle, but due to the symmetric combinations of data, 
the two responses hi {I = A, E) to the low frequency gravitational waves can be essentially regarded as that of two 
90°-interferometers rotated by 45° from each other, as shown in figure 1. We take the reference coordinate system 
(Xd, Yd, Zjj). The Z^-axis is normal to the detector plane. When we define two units vectors I and n, the response 
of the A mode is proportional to lHl~nHn and the E mode to 2n ■ H ■ I with some tensor H made from 
gravitational waves. It is important to note that there is one degree of freedom to the data combinations (A, E) j^. 
We can make the following linear combinations {dyi{(f)f),dE{4>f)) made from the original modes {dA,dE) with a two 
dimensional rotation matrix R{(j)) as 

These new modes E{(j)f)) are equivalent to the responses of two 90°-interferometers that are obtained by rota- 

tion of the original modes {A, E) with angle (j)f around the Z^i-axis. This reflects the spin-2 nature of a gravitational 
wave. We can easily confirm the following relations for the detector noise spectrum 

S'lA{<t^f)=S'iE{<^f)^S^j, = S^E. S^Eih) = 0- (2-5) 

Due to the simple relation H2.4I) this freedom of the effective rotation does not affect our final results as we see later. 
Therefore we simply use the original modes {A, E) for most of our analysis. 

Now we move to the binary confusion noise that is made by a large number of unresolved Galactic binaries. At 
the low frequency regime the response of the mode A to a monochromatic binary is given by two polarization modes 
(+, x) of gravitational waves 



resp 
III 



hA^a+FX+a'F^. (2.6) 

With the quadrupolc formula the amplitudes and are expressed as 

a+ ^K{l + cos^e,), a''^2Kcose„ (2.7) 

with the inclination angle 9i, and the coefficient K is given in terms of the distance to the binary r, the chirp mass 
rric, and the gravitational wave frequency / as 

^5/3™5/3 

K = 2 ^ 7- inff/^ (2.8) 

Two functions F^ and F^ are determined by the sky position {6, (p) and the polarization angle ip of the binary in the 
{Xd, Yd, Zo) coordinate as 

F| = i (1 + 008^6') cos 20 cos 2-0 -cos 6* sin 2(/> sin 2^;, (2.9) 

F^ = i(l + cos^6')cos2(?!)sin2?/>-|-cos6'sin2(/)cos2V'. (2.10) 

The response He to the E mode is given by a similar expression as ea. H2.6|l . We replace the azimuthal angle (j) to 
(j) — 7r/4 in eas. H2.9|l and (|2.1U|I . As in ea. (|2.2|) defined for the detector noise, we denote the Galactic binary noise 
spectrum matrix Sfj by the responses {h^, h^) of the modes {A, E) to all the unresolved Galactic binaries. With using 
ea. (|2.ti|) we sum up all the unresolved binary within Galaxy. For simplicity we assume that the binary mass function 
does not depend on the frequency / or the position in the Galaxy. This would be a fairly reasonable assumption. 
Considering the nature of the confusion noise we average over the inclination and polarization angles 6*^, -0, the spatial 
position r, and obtain the following expression 

Sfj - P / drp{r)r-^Fij{e, 0), (2.11) 
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where P is a normalization factor that is not important for our analysis. The function p{r) is the density profile of 
the Galactic binaries, and the functions Fu are defined as 



Faa = 



1 



'EE 



l+COS^e X , ^2 

cos 20 1 + (cos 9 sin 20) 



— — sin 20 ) + (cos 6* cos 2(^^^ 



Fae = Pea = ^ 



1 + cos^ 



— cos 



COS 20 sin 20. 



Using a spherical coordinate we can express ea. H2.11|l as 

Sfj = I dnB{d,cj))Fij{ 
where the anisotropy of the gravitational wave background B{t 



is defined as 



B(9,, 



P j drp{r,e, 



(2.12) 
(2.13) 
(2.14) 

(2.15) 
(2.16) 



If the source distribution p{x) for gravitational waves is same as that of the light and the interstellar absorption is 
negligible, the function B{9, 0) is also the luminosity distribution of the Milky Way on the sky. 

In reality there would be some fluctuations around the above expression (|2.11(l due to the finite number of the 
binaries. But this effect is beyond the scope of this paper. Some binaries would be close enough to have significant 
SNR above the confusion noise level, and could be removed from the background. But a rough estimation indicates 
that the typical distance to such binaries is smaller than ~ lOOpc at the low frequency region / < ImHz, and their 
subtraction would change our results only slightly 0| . 

For an isotropic binary distribution we can replace the function Fu by 



Faa = Fee 



Fae = 0, 



(2.17) 



and the noise matrix Sfj is diagonal (proportional to the unit matrix). For reference we define the noise amplitude 
^Eto that is obtained with using Fjj — 1/5 in ea. H2.11|l . This monopole spectrum Sf^^ is what has been used in most 
of previous studies. Differences between Sfj and Sj^^ due to the proper angular dependence of Fjj are the main issue 
of this paper. 

We can calculate the noise matrix Sfj{(t)f) for the new data combination (yl(0/), E{(j)f)) generated by the simple 
relation H2.4|l . It is given by the original one Sfj as 



Sfji^Pf) = R{24>f)SfjR{- 



-20/). 



(2.18) 



Let us consider an analysis for some strong gravitational wave signal {hA.x , ^e.x) from a source X, e.g. the ring 
down waveform from a merged MBH binary. When the detector noise Sfj is much smaller than the Galactic confusion 
noise Sfj, the SNR is given as 



SNR^ cx {hA,hE)*{S'' 



IlA 

He 



(2.19) 



From eas. (|2.4|l and (|2.19l) we can easily confirm that the freedom of the rotation angle 0/ is completely killed in the 
above relation for SNR. This is the reason we can forget it from the beginning. But it is instructive to take a rotation 
angle 0j = 0o so that the noise matrix Sfj{(j)o) is diagonalized as 



Then we have the following simple relation 



SNR"^ 



cx 



Ai 
As 



\hAiM\^ , IhEiMl' 



Ai 



A2 



(2.20) 



(2.21) 
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Next we evaluate the effective noise level for the two modes by averaging out the direction and orientation the source 
X. Apparently the averages for \hAx{4'f)\'^ ^-iid \hEx{4'f)\'^ do not depend on the choice of the angle and we have 
(|/iAx(</'o)P) — {\hEx{4'o)\^) ■ Thus it is reasonable to define the effective noise level 5'^^ by the following equation 

The coefficient 2 conies from the number of the modes. The inverse of this quantity is proportional to the averaged 
SNR^ for sources at a fixed distance. For ground-base detectors, such as, LIGO II, the effective distance de to a 
source that can be detected above a given SNR threshold is often used to characterize the noise level. When the 
cosmological effects for the event rate is small {e.g. ^ SOOMpc), this is a intuitive measure for a detector. But for 
the cosmological sources the situation is not so simple. Therefore we do not pursuit such a kind of measure, but use 
the effective noise Sfjj. This approach greatly simplifies the treatment of the polarization angle. 

Finally we make a brief mathematical comment on the effective noise level Sf^^. A rigid analysis is given in the 
Appendix. The functions Fij{9,(f>) are written by spherical harmonic functions YimC^^, 0) with I = 0, 2 and 4 0,0. 
This means that the angular dependence of the Galactic confusion noise felt by the detectors {A, E) is limited to 
these three modes at low frequency region. Therefore the effective noise for arbitrary detector plane is determined 
by 1 + 5 + 9 = 15 real parameters that characterize anisotropics of the binary noise (total of 5+9) and the monopole 
mode (total of 1). The monopole mode is nothing but the isotropic component S. 



III. ANNUAL MODULATION OF THE EFFECTIVE NOISE 

Due to the annual rotation of the detector plane of LISA, the effective noise Sfj^ changes with time. To denote 
the orbital phase of LISA we use the time T measured from the autumn equinox point in units of year. The actual 
date for T = corresponds to ^ 20 days after the autumn equinox day, as LISA trails 20° behind the Earth. With 
eq. (|2.11l) we can calculate the annual modulation of the noise S^j ^ at given moment T . The detector plane inclines to 
the ecliptic plane by 60°. We assume that the vertex of the cone made by the envelope of the detector planes exists 
at south as shown in figure 2. Hereafter we implicitly assume this configuration. When we take it to the opposite 
direction, the time dependence of the noise Sfjj changes by T = G.5yr from our results. In figure 2 we also define 
the ecliptic coordinate that is very useful for describing the motion of LISA. In this coordinate the directions of the 
Galactic center, the north and south Galactic poles are (6*, 0) = (-95.5°, -93.2°), (60.4°, 180.0°) and (119.6°, 0.0°) 

m 

A. Galactic model 

First we calculate the all sky r nap of the gravitational wave luminosity 5(0,0). We mainly use the Galactic 
stellar distribution model given in [ig- This model (hereafter BGS model) contains both the triaxial bulge and disk 
components obtained by fitting the near infrared COBE/DIRBE surface brightness map. In addition to these two 
components we have also studied contribution of the halo stars. But its effect is almost negligible for reasonable model 
parameters, and we do not discuss the halo component in this paper (see also [Sj)- The explicit form of the BGS 
model is given as 



p{x) = C{h{x) + f,{x)), (3.1) 



where C is a normalization constant and 



h = /oTr- — 7 — US'' (^-^^ 
(1 + a/aoY-^ 

/ -\z\/zo 

fd = [- +a- ]Rse-^"^% (3.3) 

9 9 

2 , y z-" 



R = (x2+y2)l/2. (3.5) 

Here we use a Galactic coordinate. Its z-axis is oriented to the north Galactic pole, and the direction of the cc-axis 
(the long axis of the triaxial bulge) is 20° different from the Sun-center line. In this model the solar system exists at 
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R = S.Okpc and z = 14pc. The parameters are given as /o = 624, a™ = 1.9kpc, gq = lOOpc, Rs = 2.5kpc, zq = 210pc, 
zi ~ 42pc, a — 0.27, rj — 0.5 and C — 0.6. We performed numerical integration of ea. H2.16l) and obtain the map 
B(9,(t)) as shown in figure 3 where the echptic coordinate is used. We can observe the bulge structure around the 
Galactic center. Note that, for the configuration of figure 2, the Z/j-axis of LISA passes pretty close to the north 
Galactic pole around the autumn equinox. 



B. noise map 



With the map B{9, 4>) of the gravitational wave anisotropics the effective noise map Sfj is obtained by convolution 
of the response functions Fjj. In figure 3 the noise map Sj^^j is given as a function of the orientation of the Z^i-axis 
(with angles {9, 0) in the ecliptic coordinate) of the detector. Two orientations {9, 0) and {tt ~ 9,(p + tt) give a same 
value due to the symmetry of the response functions Fjj. The Z^j-axis of LISA moves along the dot-dashed line. In 
the map the maximum value of S'^y/S'^^^ is 1.69 around the direction of the Galactic center, and minimum value 
is 0.38 around the Galactic poles. As mentioned before the Zn-axis of LISA passes pretty close to the poles at the 
autumn equinox. This means that the effective noise at the autumn equinox is almost the smallest value in the map. 
We can also see that the effective noise is higher in summer than winter. The separation between the Zo-axis and 
the Galactic center on the sky is smaller for the former. The map traces the stellar distribution along the Galactic 
disk (Milky Way on the sky), but is not as significant as figure 3. 

For the BGS model we also analyzed the deformation of the noise matrix Sfj from a following simple shape that 
is proportional to the unit matrix 



A 
A 



(3.6) 



This deviation is caused by the hexadecupole {I = 4) mode of the background B{9, 0) (see Appendix). We numerically 
calculate the ratio |Ai — A2|/(Ai -I-A2) for all the sky orientation of the Z^-axis, and found that its maximum value is 
0.152. Therefore it is not a bad approximation to put Ai = A2 = A, namely S'_f^ = = A and = = 0. In 
this approximation the noise matrix is characterized by a single quantity A(0, 0) and Figure 3 can be roughly regarded 
as its map. 

We have adopted the BGS model for the spatial distribution of the Galactic binaries. But the radial scale length 
of the disk Rg has some uncertainties • Furthermore the distribution of the gravitational wave sources might be 
different from that of the near infrared light. Here we study how the effective noise changes with the model parameter. 
We use the following exponential disk model for the Galactic density profile p{r) 

p{R,z) = noexp[-(i?/i?,)]sech2[(z/z,)], (3.7) 

with the Galactic cylindrical coordinate {R, z) and a normalization constant uq. The radial scale length Rg has 
typically Rs = 1.7 ~ 3kpc. In this model we fix the disk scale height at Zg — 200pc, and put the solar system at 
R — 8.5kpc and z = 30pc. 

In figure 4 we show our numerical results for the normalized noise level Sfj^/S^gi^ as a function of the time T . We 
choose the scale length Rs at Rs = 1.0, 1.7 and 2.5kpc. The results for Rs = 2.5kpc and the BGS model are almost 
the same. As in figure 2 the Galactic center is nearly on the ecliptic plane and also nearly normal to the two equinox 
directions. Let us first consider a simple situation that (i) the density profile is a delta function 5'^ (r — r gc) around the 
Galactic center rcc, and (ii) the center is on the detector plane with angles {9, 0) = (7r/2, 0gc) in the [Xd^Yd, Zd) 
coordinate system in figure 1. Then the mode A((j)Gc) is completely free from the Galactic binary noise and one of 
the eigen values Ai vanishes. Thus the effective noise becomes Sf^j = from ea. (|2.22l) . In reality the density profile 
p(r) has a three dimensional structure and the eigen values A^ take finite values. But the density profile are highly 
concentrated around the center rcc, and the minimum values of S^^j/S^so becomes smaller as we decrease the scale 
length Rs as shown in figure 5. The deformation |Ai — A2|/(Ai + A2) is larger for a smaller length Rs- 

As the effective noise level Sfj^ changes with time T, it would be useful to define a time average of the noise. The 
most natural definition from the standpoint of the signal analysis would be the following one 



5ff/-(T^/ ■ (3-8) 



This quantity is inversely proportional to the all sky average of SNR^ for sources at a fixed distance as the effective 
noise Sf^^. We have Sfj^/S^so — 1-54 and 1.15 for Rs = 1.0 and 2.5kpc respectively. Therefore the traditional 
estimation based on the monopole mode gives a fairly good result for the time averaged noise. 
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IV. SNR AS A FUNCTION OF THE SKY POSITION 

So far we have discussed how the binary noise changes with the orientation of the detector. In this section we 
first analyze the dependence of SNR on the sky position of a monochromatic source with a fixed distance and lyr 
integration time. We denote the source direction (0s,<^s) in the ecHptic coordinate. With the rotation of LISA, the 
noise matrix Sfj{t) and the response function F^"^ both change with time. At a given moment, the direction of 
the binary (0_d,(/)_d) in the rotating detector coordinate (X]j,Y]:), Zjj) can be formally expressed as {dD{Os, 'Ps,t), 
(pniOs, 4's,t)). Then the quantity SNR^ with averaging over the source orientation (polarization and inclination) is 
proportional to the following time integral 

G{e,,cf>s)^ I Y.^"i^D{es,<l>s,t),M9s,ci)s,t)){sfj{t))-^dt (4.1) 

J lyr J J 

Actually the polarization angle i/j in eas. H2.9|l and H2.10|l also changes with time. But its average commutes with the 
time integration. 

To deal with the complicated time dependence of the response function Fjj and the noise matrix Sfjit) we use the 
following three prescriptions with 

(1) the monopole noise model Sfj{t) — Sf^^Sij and the angular averaged response functions Fjj — Sjj/b, 

(2) the monopole noise model Sfj{t) — S^^^Sij but including the time dependence of the response functions Fjj, 

(3) including the time dependence for both the noise matrix Sfj{t) and the response functions Fij. 

We denote the integral G with above prescriptions as d, G2 and G3 respectively. The method (2) has been 
sometimes used for study of the signal analysis that would be valid at higher frequency regime where the detector 
noise dominant j3| • As the configuration of LISA is symmetric around the ZE-asds of the ecliptic coordinate (figure 
2) and the freedom of the (/)/-rotation is irrelevant here, the result 6*2(^8, ips) should depend only on the zenith angle 
9s- We can derive a simple analytic expression for the ratio G2/G1 as follows 

G2(6',.) _ 15605- 2300 cos 26*^ - 185cos46is 
Gi ^ 16384 ■ ^ ■ ' 

The function G2{0s)/Gi becomes maximum 1.08 at 0s = 7r/2 and minimum 0.80 at 9s = and tt. Due to the angular 
average of the response function by the annual rotation of LISA, the simple prescription (1) gives fairly good results 
for the isotropic noise. We can easily confirm that the angular average of the ratio G2{9s)/Gi becomes unity. 

Next we calculate the ratio G3{9s, (t>s)IGi by numerical integration of the expression 1)4. l|l . The results for all sky 
direction are shown in figure 6. We used the same Galactic model as figure 3. The ratio becomes maximum value 
1.56 around the Galactic poles and the minimum 0.86 (not at the direction of the Galactic center). The ratio G3/G1 
is a convolution of the sky dependence of the response function Fij and the time varying noise Sfj{t) whose profile is 
shown in figures 4 and 5. To interpret figure 6, let us regard the expression (|4.1|l as a integral of the function Sfj{t)~^ 
with the weight Fij. The function Sfj{t)~^ becomes maximum when the Z^i-axis of the detector is nearly oriented 
to the Galactic poles as shown in figures 3 and 4. At this epoch the weight function F/j is largest to the direction of 
the Galactic poles and and smallest to the Galactic disk, if we assume the matrix Sfj{t)~^ in the form (3.6) .In this 
manner the map traces the Galactic structure well. We should also notice that the overall 9s dependence of the ratio 
G3/G1 is similar to the ratio G2/G1, as easily expected. 

We can make a similar argument for a short-lived signal such as the ring down wave from a MBH binary. As 
mentioned above, the best time is the autumn equinox and the best source directions at that time are the Galactic 
poles (both Fij and Sfj{t)^^ are maximum). The worst season of the effective noise is summer when the orientation 
of the Z/j-axis is nearly on the Galactic disk. At this configuration the response function Fij takes the smallest value 
to the Galactic poles that are now normal to the Z^-axis. Thus both of two elements Fij and {Sfj)~^ are nearly 
minimized. As a results the SNR (not SNB?) at the direction of Galactic poles changes by a factor of 4 in three 
month. A factor of 2 comes from the changes of the noise and another factor of 2 comes from that of the response 
function. Both of two effects go to the same direction to amplify the time modulation. 

V. SUMMARY 

In this paper we have studied the anisotropics of the Galactic confusion noise background and its effect on the data 
analysis of LISA at the low frequency regime / < 3mIIz. In contrast to the traditional monopole approximation the 
anisotropics induce the correlation of confusion noises between the A and E modes of LISA, and the effective noise 
level depends strongly on the time. In the foUowings we briefly summarize our results. 
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The effective noise level of LISA becomes smallest around the autumn equinox, and largest in summer for the 

configuration in figure 2. The difference of the noise \J Sf^^ at these two epochs is a factor of 2. The detector plane 

of LISA at the autumn equinox is almost normal to the Galactic poles and the effective noise level is very close to 
the global minimum of the noise map given for the all possible orientations of the detector. The time average of the 
effective noise is different from the traditional monopole approximation by less than 50%. 

We have also analyzed the dependence of the SNR on the sky position of sources with a fixed distance and averaged 
orientation. For a long lived source the dependence after lyr integration is weak, and fluctuations on the sky are 
~ 30% level. This is mainly due to the angular average of the response function by annual rotation of LISA. For a 
short lived source the SNR at a given direction changes strongly with time. For example the SNR of a source at the 
direction of the Galactic pole changes by a factor of '--^ 4 in three months. A factor of 2 comes from the modulation 
of the confusion noise level and another factor of 2 is from that of the angular response function of the detector. 
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APPENDIX A: SPHERICAL HARMONIC EXPANSION 



The noise matrix Sfj is determined by the orientation of the coordinate system Kd] {Zd^Yo^Zo) (as defined in 
figure 1) that would change with time. In this situation the spherical harmonic expansion for a fixed coordinate 
system Kq\ (Xq, Fq: ^o) (e-ff- the ecliptic coordinate in figure 2) would be useful 0- We can related two systems Ko 
and 7^0 by the Euler angles (a, /?, 7). Here (a, (i) is the direction of the Z^i-axis in the Kq system, and 7 is essentially 
corresponds to the freedom of the (/)/-rotation in ea. H2.4|l around the Z/j-axis. We denote the noise matrix in the 
Kd system as Sfj{a,f3,'-f). In the configuration of figure 2 we can parameterize the motion of LISA with a = 27rT, 
/3 = — 7r/3 and 7 = —2t:T + 70 with some constant 70 and the orbital time T. 

The angular dependence of the gravitational wave intensity fea. (|2.16f) ') i3(0o,^o) is decomposed by the spherical 
harmonics Yi„i{9a, (l)o) as 

B{eo,M = Blr,rYirn{0O, M ■ (Al) 

Here the coefficients Bim is given as 

Bi„, = {lm\B) , (A2) 

where we have used the traditional notation for the inner product (a|6) = J d£La*{d{)^ (j)o)b{9o, (j)o), and an abbreviation 
\lm >= \Yim{Oo, (po) >■ Our goal is to write down the matrix Sfj{a,P,j) with using the coefficients Bim- 
First the matrix Sfj{a,P,'j) is formally given as follows 

Sfjia, P, 7) = {Fij{eo, M\U{a, /3, 7)"'|i?(0o, 0o)) , (A3) 

where U{a,(3,"f) is the rotation operator with the Eular angles a, (3 and 7. Using the identity X^im 1^*^ ^'^1 ~ ^ 
for the complete sets \lm >, we have 



SPj{a,P,j) - ^^(^^/j(^?o,0o)|?m)(/m|[/(a,/3,7)~'|?W)(rm'|B(0o,0o)), (A4) 

Im I'm' 

= (^"(^o,M\lm){lMU{a,P,j)-^\lm'){lm'\B{eo,M)- (A5) 

Imm' 

Next we calculate the coefficients Cijim = {Fij\lm) for the angular dependence of the response functions. We have 
Cijim = for / 7^ 0, 2 and 4. The foUowings are all of the non-vanishing coefficients; 

Caaoo = Ceeoq = —7—, (A6) 
5 

4 /tt 

CaA20 — CeE20 — -jJ-^, ( A7) 
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Caa4o — Ceeao — 77:7, (A8) 
105 

1 /IT 

CaA44 = CaA4-4 = —CeE44 ^ —CeE4-4 = -^J -jQ, (A9) 
1 / TT 

CaE44 = -CaE4-4 = g Y (AlO) 

Then ea. ljA5|l is simplified to 

Sfjia, P, 7) = Boo (00|f/(a, (3, 7)"'|00) C/joo + ^ (20|C/(a, /?, 7)"'|2m) C/J20 

+ J2 (40|C/(a, /3, 7)"' |4m) C/J40 + ^ S4,„ (44|C/(a, /?, 7)"' |4m) C/J44 (AH) 

m m 

+ (4 - 4|C/(a, 7)"' |4m) C/J4-4. 

m 

For further calculation we use the relation between the matrix elements (^lm\U{a, l3,'^)^^\m''j and the spin-weight 
spherical harmonics sYim{a,P) as 



(A12) 



After some calculation we finally obtain the explicit forms for the matrix ^^^(q;, f3, 7) as 



Saa = ^Boo + ^B2mY2m{l3, o) + ^BimYim{P, o) + ^^B^^ [-4Y4mif3, a)e*'^ + 4Y4m{l3, a)e-^'^'lAl?,) 
o 00 oio 9v70 



o 00 oiiD 9\/'0 



SIe = -^^4™ [_4r4„(/3,a)e4^T-4Ki™(/3,a)e-4'^] , (A15) 
9v 70 



where the summation with respect to the index m is implicitly assumed. These symmetric expressions would be quite 
useful to study how we can estimate the anisotropics Bim form the data Sfj. The correlation analysis can be 
used to extract the hexadecupole (/ — 4) mode and the freedom of adjusting 7 would be important here. Detailed 
studies on this issue would be presented elsewhere. 

We can easily confirm that both tr5'|j and det Sfj do not depend on the angle 7, as expected. We can derive some 
expressions for the eigen values Ai and A2 of the noise matrix Sfj. With notation {a{[3,a))av for angular average of 
a function a(/3, a) over a unit sphere we have 

^(Al+A2)a. = ^^00, (A16) 

2 5 

J((A4 + A2)^)a. = ^^Blo + ^ E + TT^b E l^4™P, (A17) 



4^^ ^ ^' 25 ^" 352 315 



^(Ai-A2)V = -|^E|B4™p. (A18) 
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FIG. 2: Configuration of LISA and tlie definition of tfie ecliptic coordinate {Xe, Ye, Ze)- The Z_B-axis is normal to the ecliptic 
plane and X_E-a:xis is oriented to the autumn equinox. The direction of the Galactic center is 9 = 95.5° and = —93.2° in 
this coordinate. With this configuration the azimuthal angle for the orientation of the^_D-axis and that for the direction of the 
detector (from the Sun) differ by 180°. 
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FIG. 3: The gravitational wave luminosity -B(S, 0) on the sky with the ecliptic coordinate. The BGS model is used for the 
Galactic binary distribution. The direction of the Galactic center (OgcPgc) is given by the triangle, and the Galactic poles 
by circles. The contours correspond to \ogiQ{B{9,(f>)/ B{9gc-,4>gc)) = —2 (white) to —3 (black) with interval 0.2. The Z_D-axis 
of LISA moves along the dash-dotted line. 
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FIG. 4: All sky effective noise map in the ecliptic coordinate. The effective noise S^jj/Sfsa is plotted as a function of the 
orientation [6, 0) of the ^_D-axis of the detector. Contours correspond to S^jf/Sf^^ = 0.4 (white) to 1.6 (black) with 0.2 interval. 
The dotted line is the ecliptic plane. The Zo-sxis of LISA moves along the dash-dotted line. The direction of the Galactic 
center and Galactic poles are given by the triangles and circles. The Galactic disk plane is given by the long dashed line. 
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FIG. 5: Annual modulation of the effective noise S^jj/S^^^ as a function of the orbital phase T of LISA. The solid line is for 
the BGS model. Other three curves correspond to the different scale length Rs of the Galactic distribution l|3.7^ . 
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FIG. 6: Directional dependence of SNR^ in the ecliptic coordinate after lyr integration. The left panel is the ratio Gsif), 4>)/Gi. 
The right small panel is G2{0)/G\ that does not depend on the azimuthal angle cj> as eq.(4.2). Contours correspond to 
S^ff/Sfso = 0.9 (black) to 1.5 (white) with 0.1 interval. The directions of the Galactic poles are given by the circles. The 
Galactic disk plane is given by the long dashed line. By definition this line is 90° apart from the poles. This map does not 
depend on the possible two configurations of LISA. 



